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1. Introduction 

The first examples of naked singularities occuring in gravitational collapse which were 
constructed with a view to testing and refining the cosmic censorship hypothesis were 
given in The singularities in question arise in the collapse of an inhomogeneous 
dust sphere. One can consider the dust sphere to be foliated by infinitely thin shells 
of matter; the singularities occur when outer shells overtake inner ones, and are known 
as shell-crossing singularities (SCS). These singularities are not considered to be serious 
counter-examples to the cosmic censorship hypothesis for resons which we review now 
and will return to below. 

The first and most convincing reason is that the fluid matter model is not 
appropriate for the study of gravitational phenomenon on the smallest scales. The 
fluid model is a macroscopic approximation which works well when it can be interpreted 
as representing the smoothed out behaviour of matter fields in a region of space-time, 
as for instance in cosmology. Applying this model all the way down to a scale where 
two infinitely thin shells of fluid intersect is not appropriate. This view of the matter 
model is now part of the cosmic censorship hypothesis, see e.g. 0; the matter model 
being used must be such that it does not develop singularities in special relativity. (It 
should be pointed out that this view is not universally supported and that some authors 
consider the shel-focussing singularity in fluid models to be a valid counter-example to 
the cosmic censorship hypothesis; see e.g. 0.) In the case of the dust matter model, 
which is a macroscopic model of collisionless matter, this point of view has received 
independent support in It was shown here that the spherically symmetric Einstein- 
Vlasov system, which involves a kinetic theoretic description of collisionless matter and 
in which dust arises as a singular limit, does not admit shell-crossing singularities. Note 
that the occurrence of shell-crossing is not entirely due to the vanishing pressure; an 
SCS can also occur when the pressure is non-zero [Tj. It has also been shown that in 
spherically symmetric perfect fluids, the global visibility of a shell-crossing singularity 
is related to an unphysical situation, namely the vanishing of the sound speed at the 
singularity 

The second reason relates to the view that certain geometric properties of shell- 
crossing singularities indicate the possibility of extending the space-time beyond the 
singularity, which would then not be considered genuine. While it is generally accepted 
that geodesic incompleteness means that space-time is singular, there does not exist a 
universally accepted definition of what constitutes a genuine space-time singularity in 
the sense of boundary constructions and levels of differentiability. However an extension 
of even very low differentiability would allow one to consider an SCS as an interior point 
of the space-time rather than a boundary point (singularity). In particular, one could 
then test Clarke's notion of generalized hyperbolicity for the SCS [Hj. For example, 
it has been claimed that the metric can be written in a form which is C° and non- 
degenerate at the singularity [Zj (but see below). Also in [Zj, it was shown that the 
SCS is gravitationally weak in the sense of Tipler: Jacobi fields carried along radial 
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null geodesies running into the singularity have finite limits at the singularity. It has 
been claimed that this is evidence that an extension through the singularity may be 
constructed. However not enough is currently known about the connection between this 
definition of the gravitational strength of a singularity and the question of extendibility 
of the space-time to make this characterization useful. 

The third reason relates to actual attempts to construct an extension. In jS], a 
dynamical extension was constructed for a (very) special case of an SCS. In [H], the 
general case was treated. Motivated by behaviour observed in the Newtonian case, the 
authors analysed the possibility of constructing an extension through the singularity 
into a region filled with three superimposed dust fiows. The analysis is incomplete, 
due to the awkward nature of the system of evolution equations obtained, but strongly 
indicates the existence of a solution to the extension problem. An extension along these 
lines for the Newtonian case has recently been given in PO]- A model of shell-crossing 
involving the collision of actual (distributional) shells of dust was studied in jTT] . The 
authors find that in order to single out an extension to the future of the collision, one 
must put in extra information by hand. 

So the question remains open: does there exist a natural extension through a shell- 
crossing singularity, and if so, is the extension unique? We address these questions 
here and find affirmative and negative answers respectively. The key is to introduce co- 
ordinates which cast the field equations in a form in which the shell-crossing singularity 
is replaced by a shock wave. We use SnG = c = 1. 



2. Existence of weak solutions 



In co-moving co-ordinates, the line element for spherical dust collapse is [12] 

ds'^ = -dt^ + —^dR^ + r'^{R, t)dQ'^, (1) 
where dQ'^ is the line element for the unit 2-sphere. The field equations yield 



r,= -^E + ^, (2) 

P =^^, 3 

the latter equation defining the density p{R, t) of the dust; m = m{R) and E = E{R) are 
preserved by the fiuid fiow. These functions are arbitrary (subject to certain energy and 
differentiability conditions) and constitute the initial data for the problem, set at t = 0. 
m is twice the Misner-Sharp mass. The evolution equation Q is readily solved; in the 
general case {E ^ 0) one obtains an implicit solution. From this, r ^i may be calculated. 
Regular initial data can lead to two types of singularity: that occurring when p diverges 
at r = (called the shell-focussing singularity), and that occurring when r = but 
r 7^ 0. This latter is the shell-crossing singularity. According to the implicit function 
theorem, the SCS can be represented locally by a function t = tsc{R) > 0; this can only 
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be calculated explicitly in the marginally bound case E = 0. The time t = tsriR) > 
of the shell-focussing singularity can always be calculated explicitly. See e.g. [7j. We say 
that a shell-crossing singularity occurs if tsciR) < tsriR) for some R> 0; this indicates 
that the shell-crossing singularity precedes the shell-focussing singularity. Note that if 
tsc and tsF are sufficiently smooth (in fact only continuous), then the SCS occurs at 
values of R forming disjoint open subsets. Furthermore, it is possible to identify these 
open subsets in terms of the initial data E{R), m{R) JHl- The co-ordinate freedom can 
be used to set r{R,0) = R, and so the initial density is /i := p{R,0) = m'/R^. For 
positive density in the matter filled region, we require m'{R) > for > 0, and so we 
see that p must diverge at the SCS. Note that if there exists a region tsc < t < tsF, 
then the density given by Q is negative in this region. One needs a new form of the 
solution which can be used beyond the SCS. 

On one level, the SCS can be viewed simply as a breakdown of the co-moving 
co-ordinate system: for fixed t, R ^ r{R,t) is no longer one-to-one. So we make the 
transformation {t,R) {t,r{R,t)). This transformation was used in [Zj. The resulting 
line element is 

ds'^ = (1 + E)-^[-{l - —)dt^ + 2udtdr + dr^] + r^dl]^ 

where u = \/ E + m/r. This appears to yield a C° and non-degenerate metric at the 
SCS. However this does not take account of the fact that m, E are functions of t and 
r. In particular, a space-time event on the shell crossing corresponds to (at least) two 
values of the co-moving co-ordinate i?, and so one cannot uniquely determine the value 
of m{R) or E{R) at the shell-crossing. Thus the line element above does not provide 
an extension through the shell crossing unless the functional dependence of m and E 
on the proper time and proper radius is given; this principal aim of this paper is to do 
this. The dependence in question is worked out by solving the field equations in the 
new co-ordinate system. These are simply the conditions that m, E are preserved by 
the fiuid fiow, and read 

m t — um^r = 0, (4) 
E^t - uE^r = 0. (5) 
We define the characteristics of this system to be solutions of 

I— (^) 

Then we see from Q and © that m, E are constant along characteristics and we can 
write down an implicit solution obtained by integrating (jHl). This solution holds until 
such time as characteristics cross; the characteristics are the fiuid fiow lines and so this 
is equivalent to shell-crossing. To extend beyond the SCS, we let f{u) = and 
rewrite the system as 

+ = (7) 

m^t — um^r = 0. (8) 
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This system is fundamentally non-conservative, although ((7j) is an inhomogeneous 
conservation law which admits the weak solution m, m G L°°(]R^) if 



for all test functions G C^(]R^). holds for differentiable solutions of (|7j) and 
extends the space of functions in which we may seek solutions to L°°(R^), where 
= (0, oo) X (0, oo). In particular, a weak solution need not be differentiable. Noting 
that 0, (jHl) share the characteristic speed —u, we must expect that discontinuities in 
u, m will occur at the same points. Then for (jHJ we will need an alternative solution 
concept to that used for the conservation law ((Zj). We use an approach pioneered by 
Volpert [Hj and developed by Le Floch In this approach, we seek functions 

M,m G BV{Q), where Q = and BV{Q) is the space of real-valued functions 
of bounded variation on Q, i.e. functions locally integrable on Q whose first order 
partial derivatives are locally finite Borel measures ^1]. Such functions have at worst 
jump discontinuities almost everywhere (a.e.) on Q. If u,v G BV{Q) and (nonlinear) 
g G C^(]R, M) then using a certain regularization g{u) of g{u), the expression g{u)dv / dx^ 
is a finite Borel measure which coincides with g{u)dv/dx'^ a.e. on f2. See jEl EI for 
details of these last two statements. This allows us to make sense of the left hand side 
of (jHl) when M,m are both discontinuous. We define a weak solution of (|7I), (jH)) to be a 
pair M, m G BV{Vt) which solve in the usual weak sense (i.e. which satisfy for 
all if) G C^{VL)) and which satisfy ni t ~ um^r = in the sense of Borel measures. The 
global existence of such solutions may be demonstrated as follows. This extends the 
standard derivation of weak solutions for conservation laws [TT)] . 

As noted above, E, m are constant along characteristics and so a classical solution 
is defined until such time as characteristics intersect. Recall that shell-crossing (i.e. 
intersection of characteristics) occurs in disjoint open subsets of the parameter space 
labelled by R, which is the characteristic label. So we may restrict attention to the 
case of one such open subset and treat others in the same way. Let us first consider the 
nature of the characteristics. 

Lemma 1 Let Uo,mQ G C"'^(R+) satisfy u'q > 0, ttiq > 0, where Uo{r) = u{r,0) > and 
^o('") = fn{r, 0) > 0. 

(i) The characteristic curves are C°° on Q and through each point ofQ, there passes 
at least one characteristic. Each characteristic reaches r = in a finite amount of 
time. 

(a) If the characteristics Ri, R2 (with Ri < R2) intersect in the region r > 0, then their 
slopes at the point of intersection satisfy 



(Hi) A given pair of characteristics intersects at most once and does so at some time 




(9) 




dt 
dr 



t > 0. 
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Proof: Part (i) is proven by noting that integrating (jUI) is equivalent to integrating the 
evolution equation (0) in the co-moving co-ordinate system. The characteristic R (i.e. 
the solution of (jHl) which satisfies r\^^Q = R) reaches r = at the shell focussing time 
tsF{R)- To prove part (ii), let Ui,U2 denote the values of u at the point of intersection 
on the characteristics Ri, R2 respectively. It suffices to show that U2 > u\. Let rs be the 
value of r at which the intersection occurs. Then < Ri < R2, since Ri is the initial 
value of r along the characteristic -Ri and the characteristics have negative slope. Then 







mo{R2] 


) - mo(-Ri) 




R2 


mo{R2] 


) mo{Ri) 



> ' ' ' ' + EoiR2) - EoiRi) 

K2 ril 

= ul{R2) - uliRi) > 

as required. The proof of part (iii) follows from the conclusion of part (ii); the 
characteristic Ri can only intersect R2 from below. □ 
We note also that the regions in which characteristics intersect form connected 
subsets of Q which may be written explicitly in terms of the characteristic label R as 

Z = {{r,t) eQ:r = r{t- R),tsc{R) < t < tsF{R), Ri < R < R2} 

where r = r{t] R) represents the characteristic with label R, and Ri,R2 are least and 
greatest values of R for which the corresponding characteristics meet others in Q. 

Let p be any point in Z°, the interior of Z. We can identify the left-most and the 
right-most characteristics which intersect at p. Let these correspond to values Rl < Rr 
of the characteristic label respectively. There exists an open neighbourhood of Rl 
and a time interval / such that the change of coordinates 7 : {R,t) hh^ [r{t;R),t) is a 
diffeomorphism from Nl x / onto an open neighbourhood Np of p. (The Jacobean of 
this transformation is singular only at the shell-crossing singularity, which corresponds 
to the boundary of Z.) We can then define a function ul '■ Np —>■ 'Rhj UL{r, t) = uo{R), 
where R is defined by {R,t) = 7~^(r, t). From the assumptions and part (i) of Lemma 
1, we see that ul is on Np. This extends to a function on the interior of Z 
and a corresponding map ur is defined in a similar fashion. Part (ii) of Lemma 1 
implies that ul,ur satisfy what we will refer to as a gap property: uiir^t) < UR{r,t) 
for (r, t) G Z°. 

Let (ro, to) be the co-ordinates of that point po on the shell-crossing singularity (i.e. 
on the boundary of Z) for which to is minimal. This is the globally earliest occurrence 
of characteristic crossing, and we assume for convenience that this defines a unique 
point. (This will involve an assumption about the initial data UQ,mQ.) The existence, 
smoothness and gap property of ul,ur imply that there exist C° piecewise curves 
r = (f){t) extending through Z° from po to r = with the property that 

-ULim.t) > <p'it) > -UR{(j){t),t), 
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for all values of t for which 4>'{t) exists. The relevance of such curves is that they provide 
a means of constructing a unique foliation of fl by characteristics. Q — Z already admits 
such a foliation. In Z and to the left of $ (the graph of 0), we use the characteristics 
which approach $ from the left (those used to define ul along $). In Z and to the right 
of $, we use the characteristics which approach $ from the right. We note that all uses 
of the words 'right' and 'left' are well-defined by the gap property and the definition of 



We are now in a position to write down weak solutions of ((Zj) and (jH)). For each 
path (j) described above, we define m^{r,0) = mo(r) and then define throughout 
f2 — $ by taking m,^ to be constant along the characteristics of the foliation constructed 
around (p. Taking rriQ G C^(]R+), Lemma 1 shows that G C^{Q — with a jump 
discontinuity across $ — po- This is sufficient to guarantee that G BV{Q). We 
do likewise to obtain G BV{Q). This procedure yields BV functions m^,E^ which 
solve the field equations Q and (0) along individual characteristics. However more is 
required in order to have a weak solution. A standard argument (see for example [TH] ) 
shows that in order for u,m to satisfy (jH)) for all tjj G C^{Q), the Rankine-Hugoniot 
condition must be satisfied: 



For any function T on Q, [T] = T^—T_ where T± denotes the right and left-hand limits 
of T{r,t) as the discontinuity is approached along t =constant. Thus when ul 7^ ur, 
i.e. anywhere in the interior of Z, we must have 



The smoothness of ul, ur implies the existence and uniqueness of a solution r = (j)^{t) 
of the initial value problem consisting of (fTT|) and the initial condition r(to) = '"o- 

It remains to verify that (jH)) is satisfied in the sense of measures. This is trivially 
the case on — where the equation is solved in smooth functions. So consider the 
measure a = m^t — um^r- For any point {r,t) on the shock, this measure evaluates to 



We have used here Volpert's definition of the regularization m of a discontinuous 
function, and applied the Rankine-Hugoniot condition after evaluating the integral. 
This concludes the proof of the existence of weak solutions of the system ([7j), (jHI). We 
summarise as follows. 

Proposition 1 Let uq^ttiq G C^(0, 00) satisfy u'q > 0, > 0. Then there exists 
u,m E BV{fl) giving a weak solution of (0) as defined above with u(r, 0) = 
uo{r),m{r,0) = mo(r) and which satisfies m{x + h,t) > m{x,t) for all h > and 

all{x,t)en. □ 

The last statement here is immediate from the construction and indicates that the 
solutions have positive distributional density. Note that the solutions are global in t, 



(10) 



m = -^M(p{t),t)) + UR{^{t),t)). 



(11) 
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but that this does not imply the absence of a future singularity. Figure One illustrates 
the result of this Proposition for the initial data Eq = and 

,„„(0=/*f)' 0<.<2.25; 

I M(2.25), X > 2.25, ^ ^ 

where 

e 



M(0 



(1 + 12^-9^2 + 2^2' 
Thus the boundary of the dust sphere is initially at x = 2.25. 



3. Uniqueness considerations 

In both the marginally and non-marginally bound cases, uniqueness fails on a 
fundamental level. To see this we focus on the marginally bound case. In this case, 
E = and there is only field equation which can be written in conservative form by 
taking x = — |r^/^; 

m^t - Vmm^^ = m^t + f,x = 

where / = — |m^/^. For this form of the equation, the Rankine-Hugoniot condition is 
[m]4> = [/]. For any smooth function M(m), the solution of the pde's for M and m agree 
in the smooth domain. However the conservative form of the equation for M will not be 
the same as that for m. This leads to a different shock speed (lllj). and hence a different 
location for the shock. For example, M = Inm leads to [lnm]0 = [—2^/rn\. Thus we do 
not have a unique weak solution for given initial data; we could only expect to obtain 
such if the field equations were fundamentally integral conservation laws which is not 
the case. 

However it is possible to obtain a uniqueness result on a different level for the 
marginally bound case. (There are no global uniqueness results for weak solutions 
of non- conservative first order systems.) We note that as with entropy solutions, the 
condition which yields uniqueness has a strong physical motivation, being equivalent 
in the present case to positivity of the energy density considered as a measure. The 
proof involves a variation of Oleinik's proof of the uniqueness of entropy solutions of 
conservation laws with convex flux functions (see chapter 16 of |l(ij). 

Proposition 2 Let rriQ G L°°(R+) satisfy mo(x + h) > mo(x) for all x > and h > 0. 
Let f{m) = — |m^/^. Then there exists a unique weak solution on of the conservation 
law 

m^t + /(m),, = (13) 
which satisfies 

m{x + h,t) > m{x, t) 
for all h > and for all {x, t) G M^, and m{x, 0) = mo{x) . 
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Proof: We refer to a solution described in this proposition as a positive density 
solution. 

(i) As regards existence, the proof of Proposition 1 is easily adapted to the 
present case. The regularity of the data can be lowered by virtue of the fact that 
the characteristics in this case are straight lines satisfying 
dx 



jm. 

dt 



(ii) For uniqueness, our aim is to prove that there is at most one positive density 
solution in L°°{VL). So let m, m be two such solutions with the same initial data. Then 
each satisfies the weak solution condition 



[mip^t + f{'m)i',x) dxdt = 
for all e C^{Q), and so we can write 

= I [{m- m)ip^t + {f{m) - f{m))ip^^] dxdt 

I 

w[ip^t + G{m{x, t),m{x, t))4'^x] dxdt, (14) 



n 



where w = m — m and 

G = [ f\\m{x,t) + {l-X)m{x,t))dX. 
Jo 

The proof proceeds by attempting to replace the term multiplying w in the integrand of 
(fT^ by an arbitrary test function (p G C^{Q). The result would then be immediate, but 
the non-smoothness of G must be dealt with carefully. As mentioned above, this proof 
follows almost identically Oleinik's proof of uniqueness of entropy solutions as given in 

m- 

(iii) For any e G (0, 1), we define the smooth functions rrf = rie*m and = ri^*m 
where rj^ is the standard mollifier and * represents convolution. Then ||m'^||oo < ll^lloo 
and rrf ^ m pointwise a.e. in f2, with similar results holding for irf. The positive 
density condition (mo non-decreasing) implies that 

m;^(x,t)>0 (15) 

in Vt. 

Next, we define = G{m'^,rff ) and write (fT^ as 

= / «;(^ t + G'i)^^) dxdt + / w{G - G^)V',x dxdt. (16) 
if7 ' ' Jn 

Now let (f) £ C*^(f2) be arbitrary, and take T > such that supp(0) C (0, oo) x (0,T). 

Let V'^ be the solution of the terminal value problem 

^'t + G'^Ij% = in (0, oo) X (0, T) 

=0 on (0,oo) X {t = T}. (17) 
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This problem may be solved by the method of characteristics and admits a unique 
smooth solution. Furthermore, since |G"^| is bounded in the compact set supp(0) C Q, 
the solution ip"" also has compact support in (0, oo) x [0,T). 
(iv) We next establish the following upper bound: 

<e^||0,,|U in supp(#). (18) 

Note that 



1 







Since m > in f2, /" = — ^m^^" is negative and bounded on any compact subset of il. 
The positive density condition (fT3jl then yields < in supp(?/'''). 
Let a = eVx- From (fTTj) . we obtain 

at + G'a,^ = (1 - G%)a + eV,a:, (19) 

and we note that a = on t = T. Let s > 0. In the compact set Kg = suppi^ip'') fl {s < 
t < T}, a attains its maximum and minimum values. Let its maximum occur at (xq, ^o)? 
we must have a(xo,to) > 0. If (xo,to) ^ dKs, with to > -s, then a(xo,to) = 0. If 
s '^tQ < T and (xq, to) is in the interior of Kg, then we must have 

a,t < 0, a^x = 
at the maximum. Then from (jl9p . at (2:0,^0) we have 

> (1 - G%)a + e*o0,x 
>a + e*V,x. 
From this we obtain 



<e*" 110,. II 00 
< e^||0,x||oo- 

Similarly, at the non-positive minimum we must have either a = or 
> 0, a^x = 0. 
In this case, (111^1 yields 

< (1 - G%)a + e*V,x 
<a + e*"0,x- 
From this we obtain 

> - e^||0,x||oo- 



Hence 



|a| < e^||0,x||oo, 
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and so we obtain the s— independent bound 

IV^:.! <e^||0,,.|U 

on supp('?/''') n {s < t < T} for any s > 0, and so we may let s — >• to complete the 
derivation of (fT8|l . 

(v) It remains to show that for any G C^{Q) and 6 > 0, 



wd) dxdt 



< 5. 



In f|TH|l . take ip to be ip", the solution of (fTTjl . Then 



w(f)dxdt = — / ^(6" — G'')ijj'',j.dxdt. 
Q Jn 

The domain of integration on the right hand side may be replaced by supp(?/'^), whereon 
flTH|) applies to bound l'?/''^.! by a constant. Then the pointwise limit G — — > yields 
the result by taking e to be sufficiently small. □ 
We note that this proof carries through for different formulations of the field 
equation (fTn|) . That is, for a diffeomorphism G : M"^ —>■ M"*" with G{0) = 0, the 
corresponding equation for M = G{m) with data corresponding to positive density 
has unique weak positive density solutions. These solutions will however disagree with 
that of Propostion 2. 



4. Discussion 



The solutions presented here show that one can extend space-time beyond a shell- 
crossing singularity in a natural way; the extension arises almost immediately when 
one steps out of the co-moving co-ordinate system. The drawback is that the extension 
is not unique. Our interpretation of this fact is that, leaving aside the question of the 
matter model, shell-crossing singularities should be included in a discussion of cosmic 
censorship. This is based on the following argument originally put forward in |17j . 
Strong cosmic censorship asserts that generically, initial data evolve to give a maximal 
globally hyperbolic space-time. It has been shown that such an evolution must be 
unique ^H]- Thus the lack of a unique evolution indicates a violation of strong cosmic 
censorship. Genericity aside, this is the situation we have found here: there are multiple 
L°° weak solutions of Einstein's equation for spherical dust for certain open subsets (c/. 
[in]) of the space of initial data. 

Having extended beyond the shell- crossing singularity, one can now investigate 
the question of the global structure of the space-time and generalized hyperbohcity 
and study geodesies impinging on the singularity. This work is in progress, utilizing 
Colombeau theory JH] as has recently been done for space-times admitting conical 
singularities [201 , [21] and for hyper surf ace singularities |22j . 

The problem of uniqueness seems hard to resolve, as it would seem to involve 
deciding on which formulation of the conservation law ((7|) is 'correct'. One possible way 
to obtain uniquely the shock surface is to consider the conservation law VqT"'' = 0. 
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This was the route taken in [221; in which the authors generahze the classic results 
of Israel [25 to Lipschitz continuous metrics and so deal with a different variety of 
gravitational shocks. In our the present case we do not have Lipschitz continuity. Indeed 
rpab measure- valued and it is far from clear how to extract jump conditions from the 
conservation law. An approach based on the recently developed geometric theory of 
nonlinear generalized functions ^] may yield the required results. 

There is an interesting analogy between shell- crossing singularities and shock waves 
in gas dynamics. In each case, the singularity disappears when a more appropriate 
matter model is used (dust replaced by the Vlasov model, inviscid fluid replaced by 
viscous). In the latter case, the shocks observed in the inviscid system reflect large 
derivatives which can arise in the more realistic model. It would be of interest to see if 
this pattern is repeated in the dust- Vlasov pair, where the existence of large derivatives 
in the latter model may be to all intents and purposes physical singularities. 
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Figure 1. The diagram shows the characteristics (sohd and dashed Hnes) 
corresponding to the mass function (|12|l . The dashed hnes are the characteristics that 
emanate from the vacuum region of the initial data surface t = 0. The sheh-crossing 
singularity corresponds to the boundary of the region Z in which characteristics cross, 
the earliest point of which is roughly at {x, t) = (0.496, 3.695). The shock, shown bold, 
evolves from this point. (In the figure, this has been constructed numerically using 
a second order Euler scheme.) The solution in co-moving co-ordinates holds outside 
Z] the solution constructed in Proposition One holds throughout Q. and so extends 
through the shell-cross. 



